Soliton generation of two-component Bose-Einstein condensates in optical lattices 
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Coupled nonlinear Schrodinger equations (CNLS) with an external elliptic function potential 
model a quasi one-dimensional interacting two-component Bose-Einstein condensate trapped in a 
standing light wave. New families of stationary solutions of the CNLS with a periodic potential are 
presented and their stability studied by direct numerical simulations. Some of these solutions allow 
reduction to Manakov system. From a physical point of view these solutions can be interpreted as 
exact Bloch states at the edge of the Brillouin zone. Some of them are stable while others are found 
to be unstable against modulations of long wavelength. The solutions which are modulationally 
unstable are shown to lead to the formation of localized ground states of the coupled BEC system. 
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I. INTRODUCTION 

Recent experiments on dilute-gas Bose-Einstein con- 
densates (EEC's) have generated great interest both 
from theoretical and experimental points of view 0. 
At ultra-low temperatures the mean-field description for 
the macroscopic BEC wave-function is constructed us- 
ing Hartree-Fock approximation and results in the Gross- 
Pitaevskii (GP) equation (jy]. The latter one reduces to 
the one-dimensional (ID) nonlinear Schrodinger (NLS) 
equation with an external potential, in particular, when 
the transverse dimensions of the condensate are on the 
order of its healing length and its longitudinal dimen- 
sion is much longer than its transverse ones (see e.g. 
^ j3j). This is termed the quasi- one dimensional (quasi- 
ID) regime of the GP equation. In this regime BECs 
remain phase-coherent, and the governing equations are 
one-dimensional. 

Several families of stationary solutions for the cubic 
NLS with an elliptic function potential have been recently 
presented in Refs A, 5] and their stability has been exam- 
ined using analytic and numerical methods [E IE S 13 • 
In the quasi-lD regime, the GP equations for two inter- 
acting BECs reduce to coupled nonlinear Schrodinger 
(CNLS) equations with an external potential flOlllll l (see 
also Sec. 2 below). When the scattering lengths of the 
two components, which characterize inter-particle inter- 
actions are close to each other, the CNLS equations re- 
duce to the Manakov system with an external potential. 



In the present paper we study the stationary two- 
component solutions of the CNLS with an external po- 
tential. Several cases of explicit solutions in terms of 
elliptic functions are analyzed and their stability proper- 
ties are studied numerically. In particular, we derive a set 
of stationary solutions with trivial and non trivial phases. 
We remark that some of the solutions presented in this 
paper were also analyzed independently in Ref. 0. In 
this work, however, all components of the CNLS were as- 
sumed to be proportional to the same elliptic function. 
We extend their results in the sense that we derive solu- 
tions of CNLS whose components are expressed through 
different elliptic functions. We also investigate the role 
played by these solutions as possible initial states from 
which localized matter waves (solitons) can be generated. 
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The paper is organized as follows. In section II we 
show how to derive one dimensional equations for cou- 
pled BEC starting from the original three dimensional 
problem using a multiple scale expansion in the small 
amplitude limit. In section III we present exact solution 
of the CNLS system with non trivial phases, while in sec- 
tion IV we analyze their limits (trigonometric and hyper- 
bolic). In Section V we derive stationary solutions with 
trivial phases for both proportional and non-proportional 
components. In Section VI we discuss the physical prop- 
erties of the trivial phase solutions and show, by means 
of direct numerical simulations, how these solutions may 
lead to the formation of localized matter waves through 
the mechanism of modulational instability. In section 
VII we will discuss the present results in comparison with 
previous work and in Section VIII the conclusions of the 
paper are briefly summarized. 
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II. BASIC EQUATIONS 

At very low temperatures, when the mean field ap- 
proximation is applicable, the evolution of two interact- 
ing BECs can be described by two coupled GP equations 
(j = l,2) (see e.g. OaiU) 



--v^ + y,(r) + 



E 

1=1,2 



(2.1) 

where atomic masses of the both components are as- 
sumed to be equal, V^(r) is an external trap potential, 
and aij are the scattering lengths of the respective atomic 
interactions (other notations are standard). In the case 
when it consists of superposition of a magnetic trap pro- 
viding cigar shape of the condensate (elongated, say, 
along the z-axis) and an optical trap inducing a lattice 
potential (which is assumed to be periodic along the x- 
axes) one has (j — 1,2) 



V,{v) - '^oj][\^x^ + y^ + z^] + U{kx), (2.2) 



U{kx) 



m 

U{k{x + L)) 



(2.3) 



Although in the last expression we have imposed equal- 
ity of the optical potential for the both components, in 
a general case one has to distinguish the linear oscillator 
frequencies, cui and ^2, when considering the two compo- 
nents corresponding to the different magnetic moments. 
For example, in the experimental settings of (13j with 
®^Rb atoms — — — \/2. This fact has natural im- 
plication on the resulting form of the effective system of 
coupled ID NLS equations. Indeed, different oscillator 
frequencies means that two components are located in 
two different parabolic potentials, and thus their effec- 
tive densities are different when the number of atoms is 
equal. As a consequence, even at approximately equal s- 
wave scattering lengths, and thus for an ~ 022, the two 
components will experience different nonlinearities (the 
latter being proportional to the atomic densities). 

Another important issue to be mentioned here is that 
a cigar-shaped BEG can be viewed as a waveguide for 
matter waves. As such it is characterized by its mode 
structure. As it is well known (c.f. with the nonlinear 
optical waveguides 0|) the intrinsic nonlinearity of a 
BEG results in the mode interaction (and thus energy 
distribution among modes). If however the nonlinearity 
is weak enough, the main state of the condensate can be 
considered as a weakly modulated ground state of the 
underline linear system. As it is clear that for a two- 
component BEG the respective small parameter is the ra- 
tio between the density energy of two-body interactions. 



(hereafter Uj 



length and Nj — J j^'jpdr is the number of atoms of the 
jth component) to the density of the recoil energy 2ma^\ • 
In other words the small parameter of the problem can 



be identified as e 



< 1 (where N 



No 



is the total number of atoms). In this situation a self- 
consistent reduction of the original 3D system H2.ll) to 
the effective ID system of the coupled equations can be 
provided by means of the multiple-scale technique. Since 
the details of such a reduction have already published 
elsewhere 0, ^3 fo'^ ^ single component BEG, here we 
only outline the main steps. 

Let us first introduce dimensionless variables 



, 1 



r' = {x',r'^) = -, 
ao 

and rewrite Eqs. (|2.1|) in the form 

zV-i = [-A' + r(r') + U{k'x') + 5n|^iP + 92^^] V'l, 

(2.4) 

#2 = [-A' + n^T{r') + U{k'x') + 5i2|V^iP + 922H2?] ^2, 

(2.5) 

where k! = koi, VL = a\/a\, 



r(r') = A2.t'- 



U{k'x') = —V{k'x'), 



and gij = AirNaij /ai. Next consideration depends on the 
magnitude of k (it is assumed that U'{x)/U{x) — 0{\). 
One can distinguish three main cases: 

(i) Hi ^ 1 . Then, the periodicity modifies the spectrum 
of the underline system introducing the effective group 
velocity dispersion. Resulting equations are just GNLS 
equations without periodic potential. This is the case 
similar to one considered in for the case of a single- 
component BEG. 

(ii) k' € (say k ~ e^). In this case the periodic poten- 
tial can be considered as smoothly varying, and somehow 
can be viewed as a limit of the case considered below. 

(iii) k' = ae where a ~ 1. This is the case when 
the potential periodicity is of the order of the effective 
length of the nonlinearity. Below we concentrate on this 
last case. 

To this end we consider two eigenvalue problems 



(-A' + Aa;'2 + r^)(^i - Eiip^, 
(-A' + n{Xx'^ + r'l))^2 - E2ip2 
whose normalized ground states are well known: 



(2.6) 



Ai/4 



^2 



7r3/4 



TT' 



3/4 



(2.7) 



niidEj = r2J-i(j + 2). 

The next steps are conventional for the multiple scale 
expansion (see e.g. 0])- Namely we introduce scaled 



is the linear oscillator variables 



e"r'^ and t„ = e"i' {n 



0, 1, 2, ...) which are considered as independent and look 
for the solution of (|2.4|l , (|2.5|) in the form 



\9ii 



IAi/2 



(2.8) 



3 



with 



i'f'^ = Qj{xi,t2)'~Pj{xQ,VQ)e 



-iEjto 



J = 1,2. (2.9) 



Here Qj{xi,t2) describes slow modulation of the back- 
ground state 12.7|l due to the nonlinearity. 

Substituting H2.8|l in (|2.4|1 . (|2.5|1 . equating all terms at 
each of the e orders, and excluding secular terms, in the 
order we obtain 

^5Qi ^ ^ Viaxi)Qi + XilQiPQi + Ml'Qu 



dt2 dx\ 



(2.10) 



where 



(2.11) 



Xi = sign(gii) j \ipi\'^dr 



912 

Iffiil 



l<Pil1^2|^dr 



23/2^3/2 ' 

1 ^ n ai2 



X2 



7^3/2 \^J7 + ly |aii| 



23/2^3/2 



and its is taken into account that ai2 =021- The system 
(|2.10() . H2.11|l is a subject of our main interest. 



III. STATIONARY SOLUTIONS WITH 
NON-TRIVIAL PHASES 

After the change of notations: t2 ^ t, xi ^ x, P = 
— X, and 61.2 — — Xi,2 the system (I2.10f) . (|2.11f) takes the 
well known form: 



respectively. We notice that nontrivial phase solutions 
imply nonzero current of the matter - it is proportional 
to \ qj{x)\'^Qjx = Cj , for each of the components - along x- 
axis, and hence seem to have no direct relation to present 
experimental setting for BECs (remember that the con- 
densate is confined to a parabolic trap). Meantime a 
system of coupled NLS equations appears to be a general 
model, havin g, f or example, applications in nonlinear op- 
tics (see e.g. [U). Bearing this in mind we consider both 
types of solutions. 

An appropriate class of periodic potentials to model 
the quasi-lD confinement produced by a standing light 
wave is given by 



V{ax) — Vosn'^{ax,k), 



(3.5) 



where sn {ax, k) denotes the Jacobian elliptic sine func- 
tion with elliptic modulus Q <k < 1. Then, substituting 
the ansatz (I3.3|l in Eqs. (|2.10() . (|2.11() and separating the 
real and imaginary part we get 

qlqixx + {hql + fiql)qt - Vosn'^iax, k)qf + ujiqf = Cf, 

(3.6) 



qlq2xx + {Pql + 62g|)g2 - V"osn^(aa;, fc)g| + uj2q2 = C|. 

(3.7) 

We seek solutions for (7^, j = 1, 2 as a quadratic function 
of sn [ax, k): 

q'^ = Ajsn'^iax, k)+Bj, j = 1, 2. (3.8) 

Inserting H3.8|) in (|3.t)|) , H3.7|) and equating the coefficients 
of equal powers of sn (ax, k) results in the following rela- 
tions among the solution parameters ujj, Cj, Aj and Bj 
and the characteristic of the optical lattice Vb, a and k: 



(62 - P)W 



A.^ 



(&i-/3)W^ 



(3.9) 



dt dx"^ 



..dQ2 , d^Q2 



+ {bi\Qi\^ + l3\Q2\^)Qi 



dt dx"^ 



-Vosn'^{ax,k)Qi = 0, (3.1) 



+ {P\Ql? +b2\Q2?)Q2 



-Vosn^[ax,k)Q2 = 0, (3.2) 



B, = -^P,A„ q = a^ApjiP, - l)a-Pjk^) (3.10) 



Ujj = (1 + fc2)a2 
W 



(3.11) 



+ ^[Pibi{b2-P)-f32f3{P-h)]~k'a'p. 



We restrict our attention to stationary solutions of 
these CNLS: 

Qj{x, t) — qj{x) eyi\3{—iujjt + iQj{x) + ikq^), (3.3) 

where j = 1,2, Ko.j are constant phases, qj and Oj{x) 
are real-valued functions and 

r dx' 

^A^)=C, / j = l,2, (3.4) 

Jo qj \x } 

where Cj, j — \,2 are constants of integration. 

Following [sl we refer to solutions in the cases Cj = 
and Cj ^ as to trivial and nontrivial phase solutions. 



where j — 1,2 and 

W^Vo-ia^k", A = xiX2-X^ = bib2-f. (3.12) 

In order that our results H3.8() - (|3.11|l are consistent 
with the parametrization H3.3|l , H3.4(l we must ensure that 
both qj{x) and &j{x) are real- valued; this means that 
Cj>0 and q'j{x) > 0. An elementary analysis shows 
that this is true provided one of the following pairs of 
conditions are satisfied: 

a) A,>0, /3, <0, j = l,2; (3.13) 

b) Aj<0, l<Pj<^, j = l,2; (3.14) 
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Although our main interest is to analyze periodic so- 
lutions note that the solutions Qj in H3.3|l are not al- 
ways periodic in x. Indeed, let us first calculate explicitly 
Qj (x) by using the well known formula, see e 



iC expiici 

.g. M- 



du 



P(aw) - p{av) 
1 



p'{av) 



2xC,{av) 



1 aiau — av) 
In 



a a{au + av) 



where p, C, a are standard Weierstrass functions. 

In the case a) we replace v by ivj, set sn.^{iavj;k) 
(3j < and 



1 



ei = -(2-fc^), 



e2 = i(2fc2-l), 



e3 = -^{l+k^), 



and rewrite the l.h.s in terms of Jacobi elliptic functions: 
du sn'^ (iav; fc)sn^(aw; k) 



sn'^ (iav; k) — sn'^ {au; k) 



-f3jx - a] 



du 



Jq S'n^{au, k) — j3j 
Skipping the details we find the explicit form of Qj [x) 

du 



Qj{x) = Cj 



Aj (sn 2 [au; k) — [3j) 

i a{ax + iavj) 

-TjX + - In — : -, 

z a (ax — laVj) 



(3.15) 



These formulae provide an explicit expression for the 
solutions Qj(x,t) with nontrivial phases; note that for 
real values of Vj &j{x) are also real. Now we can find 
the conditions under which Qj{x, t) are periodic. Indeed, 
from (|3.15|) we can calculate the quantities Tj satisfying: 



which corresponds to the case (|3.13|) . Then the solution 
is given by: 

tui = (/?! - (32)biAi + (1 + k^)a^ - a^epi,{^.19) 
W2 = (/?! - P2)PA2 + (1 + fc2)a2 - a^k^(32, 
C]^a^A^^(3j{(3j~l){l~p,k^), 
Bj=~PjA„ J = 1,2. 

The second particular case is the Manakov system; it 
corresponds to the choice bi — b2 — b. The result is 



= + /32A2) + (1 + k^)a^ - a^k^Pj, (3.20) 

^a'Ap,{P,-l){l-(3,k'), 
B, =. -f3,A,, Vo = b{Ai + A2) + 2fc2a^ j = 1, 2. 

We remark that in the two-component CNLS eqs. 
(|3.6|) . H3.7I) the constants 61, 62 and (3 are assumed to 
be positive. However in our considerations we do not 
need this restrictions and our formulae are valid also for 
negative values of 61, 62 and /3. 



IV. LIMITS OF THE NON-TRIVIAL PHASE 
SOLUTIONS 

A. The limit fc ^ 1 

In this limit the elliptic functions reduce to hyperbolic 
functions. Specifically, sn (x, k) = tanh(a;). Hence in this 
limit and for repulsive BECs the solutions have the form 



qj = Aj (tanh^ (ax) - /3j ) , j = 1 , 2 , 



Cj = -a^Ap,{l-pjf, 

uj, ^ PibiAi + P2/3A2 + (2 - Pj)a^ 



(4.1) 



(3.16) 



Then Qj {x, t) will be periodic in x with periods Tj — nijUj 
if there exist pairs of integers rrij, pj, such that: 

TTl ■ 1 

-^ = -7T [av.Cioj] + ujTj]-' , J = 1, 2. (3.17) 
Pj 

where uj (and w') are the half-periods of the Weierstrass 
functions. 

Of course the trivial phase solutions considered in the 
next sections are always periodic functions of x. 

We will list also solutions for two particular choices of 
bi , &2 and P which can be viewed as singular limits of the 
generic case considered above. The first one is 



= bib2, bi ^ 62 



(3.18) 



The potential has only a single well or a single peak 
V{x) = — Vb tanh^(aa;). The consistency condition (|3.13|) 
in this case takes the form: 



a) Aj>0, /3,<0, j = l,2; 



(4.2) 



while the second one (|3.14|) degenerates and dissappears. 

The same limit combined with the condition = 6162 
leads to: 

qj{x) = (tanh^ (ax) — Pj), 

c^i = (/3i-/32)&iAi + (2-/3i)a2, (4.3) 

c^2 = (/3i-/32)/3v42 + (2-/32)a', 
Cj = -a'Ap,{l-P,)\ 

A2 = -^Ai, 1/0 = 2a2, 
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and for the Manakov case bi = b2 = [3 = b we have: 



V. TRIVIAL PHASE SOLUTIONS. 



= + P2A2) + (2 - /3j>2, (4.4) 



Vo = b{Ai + A2) + 2a\ 
Then the nontrivial phases are equal to {(3j < 0): 

/ — TT I tanh ax \ 

Ojix) = ay'—pjX + arctanh 



(4.5) 



B. The trigonometric limit 

In the hmit fc — > 0, the eUiptic functions reduce to 
trigonometric functions and V{x) = — Vq sin^(aa:;). Then 



=Aj{sin^{ax)-Pj), J = 1,2, 
= + f3ibiAi + 13213 A2, 
UJ2 = + PiPAi +(32b2A2, 
C]^o?A]i3,{ii,-\), 

b2-P , 61-/3 



(4.6) 



Ai - 



6162 - 



A2 



6162-/32' 



i.e., Vo = 61^1+62^2. The consistency conditions (13.13(1 
and H3.14|l then take the form: 



a) > 0, I3j<0, j = l,2; 

b) A, <0, /3, >1, J = 1,2; 

If we assume /J^ — 6162 then 

= Aj{sm\ax)- Pj), J = 1,2, 
c^, =«' + (- l)-''+'(/3i-/32)6iAi, 



(4.7) 
(4.8) 



C]^a'A'^P,{P,-l), 

A2^~^A,, Vo^O, 
P 



(4.9) 



and in the Manakov case 61 = 62 = /3 = 6 we have: 
q^^ ^ Aj{sin^{ax)- (3j), J = 1,2, 



ojj ^ + b{f3iAi+ (32A2), 

C^=a^A]p,{f3,-l), 

K,=6(Ai+A2). 



(4.10) 



Therefore the phase integral 1(3. 4|l equals {(3j < 0): 

9,- = - arctan I w - — tan(Q!a;) ) . (4.11) 



In this section we consider solutions of ((3.1(1 , ((3.2(1 with 
trivial phase, i.e. Ci = C2 = 0: 



,{x,t) 



-iLjjt-\-iKQj 



qj{x), 



J = 1,2, 



(5.1) 



and we will look for different possible choices for the func- 
tions qi{x) and q2{x). This type of solutions are more 
flexible and in certain cases survive reductions of the 
constants /J^ = 6162 or the limit to the Manakov case: 
61 = 62 = p. They arejilso relevant for processes in BEC 
and nonlinear optics 14]. 

In the following we shall consider the qi{x) to be ex- 
pressed in terms of Jacobi elliptic functions, i.e. we 
assume the following ansatz: qi{x) — ^iJi{x), with 
Ji{x),i = 1,2 being one of the Jacobi elliptic function 
sn [ax, fc), cn {ax, k) or dn {ax, k) and 7^ specifying both 
the real amplitudes and the constant phases in ((5.1(1 . 
Note that the CNLS ((3.1(1 . ((3.2(1 possesses the gauge in- 
variance Qj — > Qje"*"*^-^ . This allows one to fix up con- 
veniently the initial phases of both Qj (x) . In most of the 
following examples we have made this choice by requir- 
ing that 7?> 0. Direct substitution of the above ansatz 
into Eqs. ((3.1(1 . ((3.2(1 provides a set of algebraic equations 
for the parameters whose solutions furnish exact ground 
states of the coupled BEC system. 

Case 1. We start with 

(71(2;) = 7isn (ax, fc), 52(2;) = 72cn (ax, fc), (5.2) 

The functions in ((5.1(1 are solutions of 1(3.1(1 provided 
the constants satisfy the relations: 



^17? - hi -W = Q, 
hi - hli -W^O, 
/37|+wi-a2(fc2 + l) = 0, 
627I -I- 0)2 — a^ = 0. 

where for convenience we have introduced 

W = Vo- 2k^a^. 



(5.3) 



(5.4) 



From this system we can determine 4 of the constants 
in terms of the others. Let us split these constants into 
two groups. The first one: 

Gi~{6i, 62, (3, W, a, fc,} 

consists of constants determining the equations and the 
potential and we assume they are fixed. The second 
group of constants 

G2 ~ {wi, U2, 7ii 72,} 

characterize the corresponding soliton solution. Next we 
solve ((5.3(1 and express the constants G2 in terms of Gi. 
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If 7^ 6162 we get the result: 



/32 - 6162 
b2{bi-l3) 
il3^ - bib2) 



W + a^{k^ + 1), 



(5.5) 



/?2 - 



7| = A^VF 

^2 /32 _ ' 



The constraints on the constants > and 7I > 
can be ensured in two ways: 



< 0, 
W >0, 



bi> P> 62, 
bi<P < 62, 



The case when = 6162 fixes up W by 
1¥ = 0, 

and then 



(5.6) 



(5.7) 



Case 3. Here 

gi(x) — 71 cn (ax, k), q2{x) = 72dn (ax, fc), (5.13) 

The functions in H5.13|l are solutions of (|3.1|) provided 
the constants satisfy the relations: 



bill + k'^hl + = 0, 
(ill + ^'^27l + = 0, 
6171^ + ^^1 + Pil - a^ = 0, 
/37? + 6272 +^2- a^k^ = 0. 

The solution of eq. (|5.14(l gives: 



2 ^2-/3 
7i - 



/32 - 6162 



F(/32-6i62) 



(5.14) 



VT, (5.15) 



7l 



-7l: 



= a2(fc2 + 1) _ 5^^2^ a;2 = - ^272- (5-8) 
Case 2. Here 

gi(x) = 7isn (ax, fc), 92(2;) = 72dn (ax, fc), (5.9) 

The functions in H5.9|l are solutions of H3.1|l provided 
the constants satisfy the relations: 



6172 - epil - = 0, 

Pll - k^b^ii -W = 0, 
Pll+u;i-a^{k^ + 1) = 0, 
627I + 1V2 ^ a^k^ = 0. 



(5.10) 



The solution of eq. (|5.1U|) gives 
/3-52 



ll 



UJ2 



-W 

m-bi) 

fc2(/32^6l62) 

&2(/3-&i) 



ll 



bi-f} 



kH(3^-bib2) 



W, 



W + a^ik^ + 1), 
W + a'^k^. 



(5.11) 



fc2(/32-&i62) 

The case /J^ = bib2 also fixes up W by H5.7|l and then 



12 — \ u '1 ' 



(5.12) 



bi 



LUl = a^(fc2 + 1) - ^71- W2 = a^ - \/6i&27i- 

The solutions in cases 1 and 2 exclude the possibility 
to have bi = 62 and 7^ > 0. The simplest way to see that 
is to check that for hi = 62 we have 7i + 7I = from 
(|5.5|) and il + k'^12 — from (|5.11|) : i.e. either il or 7I 
must be negative. 



,2 2 , &2(/3-6i) , /?(/?- &2) 

U)2 — K OL 



W + 



fc2(/32-6i62) /32-6162 



W^, (5.17) 



This solution allows the possibility to have bi — b2 — b. 
It reduces to 



li 



W 
T+b' 



w 



72 



k^iP + b)' 



(5.18) 



-i = «^-7?(6+^), c.2 = a^P-7?(^+^ 

Obviously to have 7^ > we need to require that 

W <0. (5.19) 

Case 4. We put: 

qi{x) = 71 dn (ax, fc), 92(2^) = 72Sn (ax, k). (5.20) 
Then the corresponding sets of parameters satisfy: 

Tf= ,'f .i-'-i^- (5.21) 



k^P^-bib2)' 



2 1 2 

uji = a k — 



/32 - 6162 ' 



bi{b2-f3)W 



k^{P^-bib2)' 



(5.22) 



UJ2 = a2(fc2 + i). 



/3(b2 - P)W 
km3^-bib2)' 



The subcase 61 = 62 = ^ is impossible since from Eq. 
(|5.21(l there follows fc272 4- ^| = 0. It is possible however 
to put l3'^ — 6162 in which case: 



72 



2 _ ^2 . /^l„,2 



-7i, 



272 7 2 
oji = a fc — bill 1 



UJ2 = a 



'{k' + l)-Vhb2ll. 



(5.23) 

(5.24) 
(5.25) 
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Case 5. Let now: 

gi(a;) = 7idn (ax, A:), q2{x) — j2Cn{ax,k). (5.26) 
Then the corresponding sets of parameters satisfy: 

2 if3-b2)W 2 iP-h)W 



/32 - 6162 



2,2 ^ /?(/3-6i)W^ ^ fei(/?-62)Ty 
= + +/fc^(/3^-M2)^ 



(5.28) 



The subcase bi = b2 ~ b simphfies further H5.27|l - 
(|5:^ to 



7i 



F(/3 + 6)' 



W 

' P + b' 



W 



J (3 + b' 



W 



(5.30) 



(5.31) 



(5.32) 



Here again it is natural to consider < 0. 

Let us finally consider three more cases in which the 
two components are proportional: 91 (x) — 7(72(2;) and 
qi{x) is one of the three functions sn{ax,k), cn{ax,k) 
or dTi{ax,k). Such an ansatz imposes on the system 
(|3.1|l . H3.2|l the compatibility condition 



7^(/3 -b2) + bi- (3 + uji-LU2=0 



(5.33) 



If is fuUfiUcd the system (EU, igSJ) reduces effec- 

tively to the one-component case, which has been already 
studied; see also Section VII below. 
Case 6. We choose: 

(7i(a;) = 7isn (ax, fc), 92(2;) = 72sn (ax, /c). (5.34) 
Then the corresponding sets of parameters satisfy: 

— U02 — oP'ik'^ + 1). 
The subcase 61 = 62 = 6 simplifies further (|5.35|) to 
W 



2 2 

7i = 72 = 



Here it is natural to consider W > Q 



tJi =W2 = a (fc +!)• (5.36) 



Case 7. Assume: 

gi(x) = 7icn (ax, fc), (72(x) = 72cn (ax, /c). (5.37) 
Then the corresponding sets of parameters satisfy: 

{P-b2)W {p-b,)W 

^1 " — fl2 — n^' T'^2=--55 — ri— > (5-38) 
p-^ - 61 62 - 01&2 

wi = 0^2 = a^ -I- 



The subcase 61 = 62 = 6 simplifies further (|5.38|) to 
W 



2 2 

7i = 72 



/3-f 5 



= W2 = a^ -f VF. (5.39) 



Unlike case 6, here it is natural to consider Vl^ < 0. 
Case 8. Let here: 

(3'i(x) = 7idn (ax, fc), (72(2^) = 72dn (ax, fc). (5.40) 
Then the corresponding sets of parameters satisfy: 

7l - ,9/^2 . . ^^ I2- ,2ia2 I, I, \' y^-'^^) 



A:2(/32 -6162)' 



A:2(/32 -6162)' 

2,2 W 

(jji = LU2 = a k + j^- 



The subcase 61 = 62 = 6 simplifies further (|5.4HI to 



2 2 

7i =72 



W 2,2 vv 



Here again it is natural to consider W < 0. 



(5.42) 



VI. MODULATIONAL INSTABILITY OF THE 
TRIVIAL PHASE SOLUTIONS AND LOCALIZED 
MATTER WAVES GENERATION 

To discuss the stability of the above solutions we shall 
adopt a physical point. To this end we remark that all 
the trivial phase solutions, are periodic functions of pe- 
riod twice the period of the lattice (recall that the period 
a of potential in Ea (|3.5|) is a = 2K{k'^)/a, where Kik"^) is 
the complete elliptic integral of the first kind). The cor- 
responding wave- number of these solutions is IC — tt /a 
which is just the boundary of the Brillouin zone of the 
uncoupled periodic linear system. Moreover, one can eas- 
ily check that each component qi{x), i — 1,2, satisfy the 
Bloch condition 



qi{x + Rn 



'9^(2;), 



la, n £ N, 
(6.1) 

i.e. the trivial phase solutions are exact nonlinear Bloch 
states (note that also in the nonlinear case the concept 
of a Bloch state is well defined by Eq. (|6.1() . Although 
nonlinearity does not compromise Bloch property (this 
being a direct consequence of the translation invariance 
of the lattice) , it can drastically influence the stability of 
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FIG. 1: Initial profile of stable cn-cn solution plotted against 
the potential profile (thick curve) . The continuous and dotted 
thin curves denote, respectively, the modulo square of qi and 
92- The parameters are fixed as: = 0.8, Vo = 1, a = 1, (3 = 
.5, bi — 1.0, 62 = 0.6. The amplitudes of the two components 
are 71 = -0.414039, 72 = 0.92582. 



200 




150 



100 t 



FIG. 2: Time evolution of the first component of the cn- 
cn solution reported in Fig0 To check stability the solution 
was slightly modulated in space with a profile of the form 
0.001 cos(.01a::). A similar plot is obtained for the second com- 
ponent in Fig. Q Parameters are fixed as in Figl. 



the states trough a modulational instability mechanism. 

The possibihty that localized states of sohton type 
can be generated from modulational instability of Bloch 
states at the edge of the Brillouin zone, was analytically 
and numerically proved for a single component BEC in 
optical lattice, in the cases of one 3], two and three 



spatial dimensions |l5|. In order to explore the same 
possibility to occur also in the present periodic two- 
components system we recourse to numerical simulations. 
To this regard we have used an operator splitting method 
using fast Fourier transform to integrate Eqs. (|3.1|l . H3.2() 
with initial condition taken as one of the exact solu- 
tions derived above modulated by a long wavelength L 
(2n/k <^ T^/L) and small amplitude sinusolidal profile. 

In Fig.l we depict the initial profiles of the two compo- 
nents cn-cn solution plotted against the potential profile, 
while in Fig. 2 we show the time evolution of the first com- 
ponent of this solution in presence of a small modulation 
of the type .001 cos(.01a;). 








10 15 20 25 
X 



FIG. 3; Time evolution of the unstable sn-sn solution (notice 
that both components of the solution are depicted at each 
time). The initial amplitudes are taken 71 = —0.41404, and 
72 ~ 0.92582. Parameters are fixed as in Fig. Q and the 
modulational initial profile is taken as in Fig|5| Notice the 
emergence of coupled soliton components out of the instabil- 
ity. 

We see that the profile remains stable for long time 
(the same is true also for the other component) indicat- 
ing that the cn-cn solution is stable against small modula- 
tions. We find that, except for this, all the other solutions 
display modulational instabilities out of which localized 
states emerge. This is clearly seen in Fig. 3 where the 
time evolution of the unstable sn-sn solution is reported. 
Notice that in contrast with Fig. 2 instability develops 
very quickly (already at time t=10), out of which two 
components bright soliton states emerge, as clearly seen 
at time t = 20 (notice that the bright soliton consists of 
two coupled solitons (one for each component) one bigger 
than the other. 

In Fig. 4 and Fig. 5 we show the time evolution of un- 
stable sn-cn solutions with different amplitude ratio of 
the sn and cn components. In Fig. 4 the stable cn com- 
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5 10 15 20 25 
X 

FIG. 4: Same as in Fig|3|but for the unstable sn-cn solution. 
The initial amplitudes are taken as 71 = —0.41404, and 72 ~ 
0.92582. Parameters are fixed as in Fig0 Notice that the cn 
component is larger and more stable. At t=100 a bright-dark 
soliton is formed in the center. 




FIG. 5: Same as in Figl^but for the unstable sn-cn solution. 
The initial amplitudes are taken as 71 — 0.92582, and 72 ~ 
—0.41404. Parameters are fixed as in Fig0 Notice that the 
unstable sn component dominates and soliton generation is 
more effective. 



ponent is larger than the unstable sn one, while in Fig. 
5 we have the opposite. We see that, although in both 
cases we have instability, the case with larger stable com- 
ponent is obviously much more stable and less effective in 
creating localized states than the other. Also notice from 
Fig. 4 that a small amplitude localized state is formed 



from the sn component at time t=100 in the middle of 
the line which seems to have a character opposite to the 
one of the other component. We see indeed that when 
matter density is higher in one component it is lower in 
the other, this suggesting a sort of bright-dark coupling. 

A better characterization of all possible states arising 
from the mixing of stable and unstable components re- 
quires a more accurate analysis. It is interesting to in- 
vestigate also solutions involving dn components since 
these, in contrast with sn and cn components, have non 
zero spatial average, i.e. they are periodic waves on top 
of a constant background. 




FIG. 6: Same as in Fig^lbut for the unstable dn-dn solution. 
Initial amplitudes are 71 = —0.46291, and 72 — 1.035099. 
Other parameters are fixed as in FiglTI 

In Fig. 6 we depict the time evolution of a dn-dn solu- 
tion from which we see that it is modulationally unstable, 
leading to the formation of brigth solitons of the same 
type observed for the sn-sn case. In Figs. 7-8 we depict 
similar evolutions for the cases dn-sn and dn-cn. Also in 
this case we observe that the mixing with the unstable sn 
component is more effective than the one with the stable 
cn component in creating localized excitations of soliton 
type (the three brigth solitons formed in Fig. 7 at time 
t « 10 remain equally spaced and well localized also on 
longer times) . By increasing the cn component of the dn- 
cn solution of Fig. 8, we also find that the time evolution 
becomes more stable, as discussed for the sn-cn case. 

This analysis shows that the exact trivial phase so- 
lutions of the previous section are very useful to create 
localized excitations of two components BEC in optical 
lattice. The fact that these solutions are Bloch states 
at the edge of the Brillouin also suggests a way to cre- 
ate them in a real experiment. One could indeed start 
from a uniform density distribution of the matter in the 
potential wells, and accelerate the lattice until the state 
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reaches the edge of the Brillouin zone where the modu- 
lational instabihty take place. 




FIG. 7: Same as in Fig|3]but for the the dn-sn solution. Initial 
amplitudes are 71 = -0.92582, and 72 = 0.46291. Other 
parameters are fixed as in FiglTI 




FIG. 8: Same as in FigElbut for the dn-cn solution. Initial 
amplitudes are 71 = 0.46291, and 72 = -0.92582. Other 
parameters are fixed as in Fig^ 



VII. DISCUSSIONS 

In this section we briefly discuss the above results 
in comparison with those of Ref |l2j | in which an n- 
component NLS-type equation with external potential, 



whose strength can be different for each component, was 
also considered. Changing somewhat notations to avoid 
confusion with ours we write it down in the form: 



1 d^^Pj 



dt 2/ij dx^ 
Vj{x) = — VbjSn^(aa;, fc) 



n. (7.2) 



1, 



In Ref. [131 three types of trivial phase solutions are 
analyzed in more details for a rather special ansatz for 
in fact it is required there that all up to a stan- 
dard phase factor are proportional to the same function 
ili{x,k). This means that the systems of n equations will 
reduce to just one equation for 'ip{x,k)] the remaining 
n — 1 equations must follow as a consequence of the first 
one and a set of constraints on the coefficients Oj/, Nj, 
fij, Uj in the notations of 12]. The same argument holds 
true also for three of our solutions, cases 6, 7 and 8 which 
we added just for the sake of completeness. 

The Hamiltonian of the n-component NLS ()7.1|l is: 



H = dx 



n 

T — 



k=l 



dx 



1 " 

- ajp\iPjf\iPp 



(7.3) 



where the integration goes over one period. Let us now 
assume that our solution is of the form: 

^l;j{x) = nj{x,t)i^{x,t), nj{x,t) = e"^'^.*+^e.(^) ^iV^ 

(7.4) 

where Nj > and Qj{x) appears only in the non-trivial 
phase case and is determined by: 



dx iV,-|i/'(a;)|^ 



(7.5) 



Inserting H7.4(l into the Hamiltonian we easily get the 
following reduced Hamiltonian: 



Hied = dx 



dip 



dx 



n ^ 



V{XM^ 

(7.6) 



V{x) — vosn'^{ax,k), 



n 

Wo = J2 ajpNjNp. 



which describes the dynamics of the effective field ip{x, t). 
The result for the trivial phase solution case leads to i?rcd 
with Af_i = 0. That is why we believe that the multi- 
component effects should be analyzed by using ansatz 
more general than H7.4|) . 
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At the same time an important result of [l^l is the 
detailed, both analytical and numerical, analysis of type 
B solutions to l|7.1|l though again using H7.4|l . 

In the present paper we considered intrinsic two- 
component solutions i.e. solutions with different ampli- 
tudes. These solutions seems to have stability property 
which are not trivial consequence of the theorems proved 
in and deserve additional studies. Besides enlarging 
the set of solutions, we have also shown the role played 
by these solutions as initial states from which localized 
matter waves (solitons) can be generated through the 
modulational instability mechanism. 

VIII. CONCLUSIONS 

In conclusion, we have considered the two-component 
CNLS with an elliptic potential as a model for trapped, 
quasi-one-dimensional two-component BECs. Classes of 
elliptic, solitary wave and trigonometric solutions have 
been presented. 

From a physical point of view the solutions discussed 
in this paper are exact nonlinear Bloch states with wave 
number at the edges of the Brillouin zone. These solu- 
tions, except for the cn—cn one, are unstable under small 
amplitude and large wavelength modulations. Two com- 
ponent matter solitons arise from these unstable solution 
via a modulational instability mechanism which resemble 



the one described in Refs. [j, ll5| for single component 
BEC in optical lattices. 

Further perspectives of finding stable periodic solu- 
tions to the 2-component problem could be linked to in- 
vestigations of finite-gap solutions of Manakov system 
given in terms of multi-dimensional 0-functions |l7j and 
[l8l | and to reduction of finite-gap solutions to elliptic 
functions j^^. Interesting classes of periodic solutions 
can be also obtained as the result of reduction of the Man- 
akov system to completely integrable two-particle system 
interacting with fourth order potential |^ and . 
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